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The formulation and implementation of an optimization method called Simultaneous 
Aerodynamic and Structural Design Optimization (SASDO) is shown as an extension of 
the Simultaneous Aerodynamic Analysis and Design Optimization (SAADO) method. It 
is extended by the inclusion of structure element sizing parameters as design variables and 
Finite Element Method (FEM) analysis responses as constraints. The method aims to 
reduce the computational expense incurred in performing shape and sizing optimization 
using state-of-the-art Computational Fluid Dynamics (CFD) flow analysis. FEM struc- 
tural analysis and sensitivity analysis tools. SASDO is applied to a simple, isolated. 3-D 
wing in inviscid flow. Results show that the method finds the same local optimum as 
a conventional optimization method with some reduction in the computational cost and 
without significant modifications to the analysis tools. 
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Nomenclature 

wing scmispan 
drag coefficient 
lift coefficient 

pitching moment coefficient 
wing root chord 
wing tip chord 
design objective function 
design constraints 
stiffness matrix 
aerodynamic loads 
free-stream Mach number 
unit normal vector 
local aerodynamic pressure 
compliance, the work done by the aerody- 
namic load to deflect the structure 
free-stream dynamic pressure 
flow-field variables (state variables) at each 
CFD mesh point 
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change in flow solver field variables due to 
better analysis convergence 
change in flow solver field variables due to 
design changes 

aerodynamic state equation residuals at 
each CFD mesh point 
norm of the residual ratio, current/initial 
surface area 

semispan wing planform area 
structural deflections (state variables) 
change in deflections due to better analysis 
convergence 

change in deflections due to design changes 
wing weight 

CFD volume mesh coordinates 
location of wing root leading edge 
chordwise location normalized by local wing 
section chord 

longitudinal location of wing tip trailing 
edge 

root section maximum camber 
free-stream angle-of-attack 
design variables 
structural element size factor 
line search parameter 

operator which indicates a change in a vari- 
able 
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S.s' convergence tolerances for function and gra- 

dients 

A adjoint variable 1 for Q 

/' adjoint variable for u 

t twist angle at wing tip, positive for leading 

edge up 

Subscripts 

d deflected shape 

j jig (undeflected) shape 

■s wing surface mesh 

r wing volume mesh 

Superscripts 

* designates updated value 

gradient with respect to design variables 

Introduction 

Simultaneous Aerodynamic Analysis and Design 
Optimization (SAADO) is a procedure that incorpo- 
rates design improvement within the iteratively solved 
(nonlinear) aerodynamic analysis so as to achieve fully 
converged flow solutions only near an optimal design. 
When SAADO is applied to a flexible wing rather 
than a rigid wing, the linear Finite Element Method 
(FEM) solution is iteratively coupled with the nonlin- 
ear Computational Fluid Dynamics (CFD) solution. 
When design variables that control structural element 
size are included, it is renamed Simultaneous Aerody- 
namic and Structural Design Optimization (SASDO). 
Overall computational efficiency is achieved because 
the many expensive iterative (nonlinear) solutions for 
non-optimal design parameters are not converged (i.e., 
obtained) at each optimization step. One can ob- 
tain the design in the equivalent of a few (rather 
than many) multiples of the computational time for 
a single, fully converged coupled aero-structural anal- 
ysis. SAADO and similar procedures for simultaneous 
analysis and design (SAND) developed by others are 
noted and discussed by Newman et al. 1 These SAND 
procedures appear best suited for applications where 
the discipline analyses involved in the design are non- 
linear and solved iteratively. Generally, convergence 
of these discipline analyses (i.e., state equations) is 
viewed as an equality constraint in an optimization 
problem. From this latter point of view, the SASDO 
method proceeds through infeasible regions of the de- 
sign space which includes not only the design variables 
3, but also the state variables Q and u. A further ad- 
vantage of SASDO is the efficient utilization of existing 
discipline analysis codes (without internal changes), 
augmented with sensitivity or gradient information, 
and yet effectively coupled more tightly than is done in 
conventional gradient -based optimization procedures, 
referred to as nested analysis and design (NAND) 
procedures. 1 A recent overview of aerodynamic shape 



optimization 2 discusses both NAND and SAND pro- 
cedures in the context of current steady aerodynamic 
optimization research. 

For single-discipline design problems, the distinc- 
tion between NAND and SAND procedures is fairly 
clear and readily seen. With respect to discipline fea- 
sibility (i.e., convergence of the generally nonlinear, 
iteratively solved state equations), these procedures 
can be viewed as accomplishing design by using only 
very well converged discipline solutions (NAND). or 
as converging a sequence of discipline solutions from 
poorly to well as the design progresses (SAND). How- 
ever, the problem formulation and solution algorithms 
may differ considerably. About twenty SAND refer- 
ences are quoted by Newman et al. 1 and Newman et 
al.; J these references discuss a variety of formulations, 
algorithms, and results for single-discipline problems 
(mostly CFD applications) in the sense of SAND as de- 
fined above. For multidisciplinary design optimization 
problems, the distinction between NAND and SAND 
is somewhat blurred because there are feasibility con- 
siderations with respect, to all the individual discipline 
state equations, as well as with respect to the multi- 
disciplinary system compatibility and constraints. A 
number of the papers in Ref. 3 discuss MDO for- 
mulations and algorithms that are called SAND-like; 
however, not all of these latter MDO procedures ap- 
pear to agree with the sense of SAND defined above 
and used herein: one that does is Hef. 4. 

The computational feasibility of SAADO for quasi 
1-D nozzle shape design based on the Euler equation 
CFD approximation was demonstrated by Hon et al. 5 
and Maui.'' Application of SAADO for turbulent tran- 
sonic airfoil shape design based on a 2-D thin-layer 
Navier- Stokes CFD approximation was demonstrated 
and reported in a later paper by Hou et al. 7 Both of 
these application results are summarized and briefly 
discussed in Ref. 1. The application of SAADO for 
rigid ICD wing design based on the Euler CFD ap- 
proximation was presented in Hef. 8. These SAADO 
procedures utilized quasi-analytical sensitivity deriva- 
tives obtained from hand-differentiated code for the 
initial quasi 1-D application 5 " and from automatically 
differentiated code for both the 2-D airfoil applica- 



tion' and the 3-D rigid wing applicatic 



Different 



optimization techniques have also been used in these 
SAADO procedures. 

The extension to multidisciplinary analysis with 
shape design variables only was presented in Ref. 9. 
Our initial results from SASDO are given in this pa- 
per. The analysis problem, the objective function, and 
the constraints are the same as those used in Ref. 9. 
That is, changes in design variables are sought to 
produce improvement in the lift-to-drag ratio of a sim- 
ple wing subject to both aerodynamic and structural 
solution-dependent constraints. These constraints are 
the difference between the lift and weight, the pitch- 
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ing moment coefficient, and the compliance, a function 
representing work done by the aerodynamic load to 
deflect the structure. There are also geometric con- 
straints. 

The flexible wing studied here is formulated as a 
static aeroelastic problem. Similar problems have 
been used as examples in Refs. 10-15 to study var- 
ious solution strategies for multidisciplinary analysis 
and optimization. In particular, Arian 1J analyzed the 
Hessian matrix of the system equations to derive the 
mathematical conditions under which the aeroelastic 
optimization problem can be solved in a "loosely" cou- 
pled manner. The multidisciplinary research of Walsh 
et al. 1314 emphasized engineering aspects of integrat- 
ing high fidelity disciplinary analysis software and dis- 
tributed computing over a network of heterogeneous 
computers. The aeroelastic analysis results of Reuther 
et al. 15 were verified with experimental data. 

Only a limited amount of literature related to aero- 
elastic problems has elaborated on the coupled sen- 
sitivity analysis. Kapania, Eldred and Barthelemy: 1 " 
Arslan and Carlson; 1, and Giunta and Sobieszczanski- 
Sobieski 18 derived global sensitivity equations (GSEs); 
some matrix coefficients in these GSEs were evaluated 
by finite differencing. Guinta 1! ' later introduced modal 
coordinates to approximate the elastic displacement 
vector in order to reduce the size of the GSE. Newman, 
Whitfield, and Anderson 211 used the complex variable 
approach to obtain the aeroelastic sensitivity deriva- 
tives, whereas Reuther et al. 1:j employed the adjoint 
variable approach to derive the aeroelastic sensitivity 
equations. A mathematical study of the coupled non- 
linear, incompressible aeroelastic analysis and sensitiv- 
ity analysis problems has been given by Ghattas and 
Li." 1 Recent results on aeroelastic sensitivity analysis 
and optimization can be found in Refs. 22-24. Partic- 
ularly, Maute et al. 23 and Hou and Satyariarayana - ' 4 
explicitly formulated the deflection update and the 
load transfer between the separate flow and structures 
solvers as part of the coupled sensitivity equations. 
In the present study, the coupled sensitivity equa- 
tions are constructed by differentiating the aeroelastic 
state equations and solving them by a Generalized 
Gauss-Seidel (GGS) method. 12 The present SASDO 
concept is very similar to that of Ghattas and others, 
Refs. 4,21,25,26, but differs in the implementation 
details as described later. 

Problem Description 

To evaluate the efficacy of the SASDO procedure for 
a problem involving multidisciplinary analysis, it, is ap- 
plied herein to a simple, isolated, flexible wing. The 
wing shape consisted of a trapezoidal planform with a 
rounded tip. It was parameterized by fifteen variables; 
five described the planform, and five each described 
the root and tip section shapes. A schematic of the 
wing and its associated shape parameters is shown in 





Fig. 1 Description of seinispan wing parameteri- 
zation. 

Fig. 1. The baseline wing section varied linearly from 
an NACA 0012 at the root to an NACA 0U08 at the tip. 
The wing structure consists of a skin, ribs, and spars. 
The ribs and spars consist of shear webs and trusses. 
Six spanwise zones of the structural mode] are defined 
as depicted in Fig. 1 . The relative sizes of the skin and 
web thicknesses and the truss cross section areas are 
fixed within each zone. Each zone is assigned a param- 
eter T which multiplies all the thicknesses and areas 
of the structural elements in that zone. The specific 
parameters selected as design variables in the sample 
optimization problems are identified in the section en- 
titled Results. The objective function to be minimized 
was the negative of the lift-to-drag ratio, —L/D. Both 
coupled solution-dependent, and geometric constraints 
were imposed. 

The solution-dependent constraints were 

• lower limit on the difference between the total lift 
and the structural weight, (";, * ^ * <7-v. — lb 

• upper limit on compliance, P — § pit ■ nds 

• upper limit on pitching moment, ('„,, in lieu of a 
trim constraint 

The purely geometric constraints were 

• minimum leading edge radius, in lieu of a manu- 
facturing requirement 



side constraints 
variables 



(bounds) on the active design 
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SASDO Procedure 

Formulation 

The flexible SASDO approach formulates the design 
optimization problem as follows: 



mill F(Q. X<i,,(rf, u). it. d) 

/i.Q.u 



subject to 



!li(Q,Xdv(-1.u),u,3)<i); 



1,2. 



where the flow field Q and the structural deflection u 
are solutions of the coupled flow equation 

R(Q.X dr (d.u),d) = U (3) 

and the finite element structural equation 

K(X J (d).d)u = L(Q.X tiv {d.u)) (4) 

The deflected volume mesh, A",;, , is determined by 
the deflected surface mesh. A",/,, as A'<<, = A,/, ( A',/., ) . 
The deflected surface mesh is a result of the jig shape 
augmented by the elastic deflection, a. as A~,t, — 
Xj(:i) + u. The two disciplines are coupled through 
the deflection, (/, and the load, L. 

Recall that Q, R, and A rf! , are very large vectors. 
This formulation treats the state variables, Q and u, 
as part of the set of independent design variables, and 
considers the state equations lo be constraints. Be- 
cause satisfaction of the equality constraints, Eqs. (3) 
and (4), is required only at the final optimum solution, 
the coupled steady-state aero-structural field equa- 
tions are not converged at every design-optimization 
iteration. The easing of that restriction is expected 
to significantly reduce the excessively large compu- 
tational cost incurred in the conventional approach. 
However, this advantage would likely be offset by the 
very large increase in the number of design variables 
and equality constraint functions, unless some reme- 
dial procedure is adopted. 

Approximations 

The SASDO method begins with a linearized design 
optimization problem solved for the most favorable 
change in the design variables, Ad, as well as for the 
changes in the state variables, \Q and Air. that is, 

>«!« F(Q,X dv (d.u),u.d) 



«Q av + V Ou + ox d „ ox d , ) ^" 



+ 



, / OF OX,,,, yf , of] a, j 
+ \oX dt . UX d ,- X J + 0,i) ^> 1 



subject to inequality constraints 






/= 1,2. 



and equality constraints 



+ 



/ JUL. Mir y, , OH \ A H - I) 
\0X d „ 0X d , A J + i)d J ^ - - U 



J) 



mid 



KUKjbi),^ - L(Q, X dv (d, «)) - f^AQ 
+ (l< 



. ox. ox,, v ux 



A.i+(^«)A/? 



oi. ox d ,, 

f)X dt . 0X d . 



(8) 



'.; 



Note that Eqs. (5) through (8) are linearized approxi- 
mations of Eqs. (1) through (4), respectively. 

In this formulation, neither the residual of the non- 
linear aerodynamic field equations, R(Q.X,d). nor 
that of the structures equation. Ku — L, is required to 
be zero (reach target) until the final optimum design is 
achieved. The linearized problem of Eqs. (5) through 
(8) is difficult to solve directly because of the number 
of design variables and equality constraint equations. 

Dmct dtffi initiation method 

One way to overcome this difficulty is by the di- 
rect differentiation method. In this method AQ. Au, 
and Eqs. (7) and (8) are removed altogether from 
the linearized problem by direct substitution. This 
is achieved by expressing XQ and Au as functions of 
Ad. 



AQ = AQi + AQ.Ad 
An = A«i + Au'iAd 



(9) 



where vectors AQi and Au\ are corrections in the 
aeroelastie solution due to the improvement of cou- 
pled aeroelastie analysis, while matrices AQ 2 and Au-j 
are corrections due to changes in the design variables. 
These vectors and matrices are solutions of the follow- 
ing coupled sets of equations, obtained from Eqs. (7) 
and (8): 



dL 



^^Av, 



SA',/„ SA',i, 



A-Ajh = j^AQ, + T< 



OL 0X C 

' ox, 



10) 



Alii 



where, for the linear FE.M, Ku - /, = at every iter- 
ation, and 

Mao, 4- 0Ii aA -"- i V 4- \ii.a i <>R — n 



KAll; 



§^ A Q? + oxv 



ox 



_ (ok. v + <±K) u 

\dX,- K J + Orl J u 



ox d 



(-Yj+Au,) 



(11) 



Note that the number of columns of matrices AQ 2 
and Ami is equal to the number of design variables. 
d. Thus the computational cost of Eq. (11) is directly 
proportional to the number of design variables. 

A new linearized problem with Ad as the only design 
variables can be obtained by substituting Eq. (9) into 
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Eqs. (5) and (6) for AQ and Au. 



mm 

A/3 






dF , 9F 9; 
9« ' 9A' rf „ 

9F 



dQ ^- 9A dl , dA,/, J 
y^Aj,, 9A.j, ' du J - dp J 



;i2) 



subject to 



+ {l^AQ 2 + ^-|A^(A'; + A« 2 ) 

+ ^Aw> + %}a;?<0; /=1.2 ttj 

(13) 
The appearance of AQi and A«i in the formulation 
indicates the difference between the SASDO (SAND) 
method and the conventional (NAND) aerodynamic 
optimization method. The AQi and Au\ not only 
constitute a change in Q and u, but also play an impor- 
tant role in defining the objective function of Eq. (12) 
and the constraint violation of Eq. (13). We can di- 
rectly solve Eqs. 10 for AQi and Au\, and, in fact, 
in previous SAADO applications that is how AQi was 
determined. However, since AQi and A«i, as shown 
in Eq. (10), represent a single Newton's iteration on 
the coupled flow and structures equations, it is possi- 
ble and less computationally expensive to approximate 
their influence on the solution Q and u by several New- 
ton's iterations of the coupled aeroelastic equations. 
That is, AQi and A«i are not determined explicitly, 
but rather the first three terms of Eqs. 12 and 13 are 
viewed as updated values of F and #,. Note that the 
terms in parentheses in Eqs. (12) and (13) are ap- 
proximated gradients of the objective and constraint 
functions. Once established, this linearized problem 
can be solved using any mathematical programming 
technique for design changes, A3- Results presented 
in this paper are computed using this direct differen- 
tiation approach. 

Adjoint method 

An alternative way to remove AQ and Au from the 
linearized problems, Eqs. (5) and (6), is the adjoint 
variable method. The adjoint variables, A and //, can 
be introduced as the solutions of 



(M) T A 



dQ 
A'// 



9L}T 
dQ' ' 

dx dl , 

9AV, 
T 



• dF \T 



( ah 

\dX A 



(14) 



9A' rfl , ifA'd 



dF flA'rf,, ) _ <±F_ 

du 



so as to rewrite the expression of the objective function 
in Eq. (5) in terms of A, //, and A3 as 



F(Q, -\'dv{3, u), u. 3) + X T R + /i T (A'« - L) 

dF dX, u . yl , 9f \ 
9A'd, dX d , - v j "•" dji ) 



(0 



j 






(15) 
Note that the terms in the brace represent the gra- 
dient of the objective function. The terms X T R and 
/< T (A'u — L) indicate the effect on the design opti- 
mization formulation due to errors in the aeroelastic 
analysis. Furthermore, in the linear sense, the adjoint 
variables and the solution errors can be related by the 
following equations: 



and 



9F _ h t 



^ = (A-«-I) T 



;i6) 



:it: 



These equations have been mentioned, for example, 
by Pierce and Giles 27 and Venditti and Darmofal 28 
for aerodynamic problems. 

In the typical optimization problem there are many 
design variables. When one can also pose the op- 
timization problem such that there are only a few 
output quantities for objectives and constraints or. in 
the extreme, combine the constraints and objective 
function into a single cost function, the adjoint ap- 
proach to sensitivity analysis has the advantage that 
the adjoint solutions are independent of the number of 
design variables. jy However, when the disciplines are 
loosely coupled, this approach is impractical since the 
coupled sensitivity analyses would require an adjoint 
for each disciplinary output being transferred, i.e., the 
discretized loads and deflections. In a tightly or im- 
plicitly coupled multidisciplinary analysis the adjoint 
approach may prove practical since this system is anal- 
ogous to a single discipline. 

Line Search 

A one-dimensional search on the step size parame- 
ter 7 is then performed in order to find the updated 
values of A3, AA', AQ, and Au. Given the search 
direction A3 determined by either the direct differen- 
tiation method (Eqs. 12 and 13) or the adjoint method 
equivalent, this line search functions to adjust its mag- 
nitude so as to simultaneously ensure better results for 
both design and analysis (converged solutions). The 
step size parameter 7 plays the role of a relaxation 
factor in the standard Newton's iteration. The search 
procedure employed solves a nonlinear optimization 
problem of the form 



mm F(Q*,X*,u*,jT) 



;i8) 
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subject to 



n(Q*.x\u\ir)<o- 



i= 1,2, 



H(Q',.\',3') =0 



h'(.\",3')u m 



(19) 



(20) 



(21) 



L(Q*,A") 

where step size -y is the only design variable. Again 
it is noted for emphasis that the equality constraints. 
Eqs. (20) and (21), are not required to be zero (reach 
target) until the final optimum design; violations of 
these equality constraints must simply be progressively 
reduced until the SASDO procedure converges. 

The updated Q' and w* can be viewed as Q* = Q + 
AQ and u* = u + An. where AQ and Am satisfy the 
first order approximations to Eqs. (20) and (21). That 
is, AQ and A« are the solutions of Eqs. (7) and (8) 
where, in Eq. (9), A3 is replaced by A3" = -)A3. 
Consequently, Q* = Q + AQ, + - ) AQ- 2 A3 and «* = 
(«+A«i +-) Au->A3 are readily available once -, is found. 
The A()i terms appearing in the above SASDO for- 
mulation are due to better convergence of the coupled 
analysis, whereas the A()-j terms are due to changes in 
the design variables. In fact, AQ? and A«j approach 
the How field and deflection sensitivities, Q' and u'. as 
the solution becomes better converged. 

Implementation 

The following pseudocode shows algorithmically 
how the method was implemented. 
set initial analysis convergence tolerance. 5 
set initial solution vectors, Q and a, 
set initial design variables. 3 
do until converged 

1. solve coupled aeroelastic analysis, Eqs. (3) &: (4), 
partially converged to = 

2. compute F and g 

3. solve coupled aeroelastic sensitivity analysis. 
Eq. (11), partially converged to s' 

4. compute A3 terms of Eqs. (5) A: (6) 

5. solve optimization problem Eqs. (5) k (()) for A3 

(5. solve Eqs. (18) through (21) for line search pa- 
rameter, ■) 

7. update 3, u, and Q 

8. tighten analysis convergence tolerance, e — e * Z, 
0< Z < 1 

enddo 

This pseudocode is similar to that used in the 
Biros and Ghattas 25 SAND approach. Specifically, 
both approaches use a Sequential Quadratic Program- 
ming (SQP) method to solve the design equations 
(step 5) and an approximate factorization method to 
solve the system equations (step 1). Step 3 above uses 
an incremental iterative method with approximate 
factorization to solve for derivatives in direct mode 
rather than as a solution of the adjoint equation of 
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Fig. 2 Diagram of SASDO procedure. 

Biros and Ghattas. 25 In addition, the line search step 
(step 6) and the convergence tightening step (step 8) 
were not included in the Biros and Gliattas method. 
A schematic of the present SASDO procedure is shown 
in Fig. 2. The dashed box, labeled "Partially Con- 
verged System Analyis," depicts the coupled analysis 
iteration loop. Steps 1 and 2 of the pseudocode; that 
labeled ••Partially Converged Sensitivity Analysis'" 
depicts the coupled derivative iteration loop. Step 3; 
that labeled "Partially Converged Design" depicts the 
design steps. Steps ">-8 of the pseudocode. Specific 
computational tools and methods used to perform 
the tasks depicted by the solid boxes in Fig. 2 are 
identified in the next section. 

Computational Tools and Models 

Major computations in this SASDO procedure are 
performed using a collection of existing codes. These 
codes are executed by a separate driver code and 
scripts that implement the SASDO procedure as just 
discussed. Each code runs independently, some simul- 
taneously, on different processors, and the required 
I/O transfers between them, also directed by the 
driver, are accomplished via data files. 

The aerodynamic flow analvsis code used for this 



study is a version of the CFL3D code 



Only Eu- 



ler analyses are performed for this work, although the 
code is capable of solving the Navier-Stokes equations 
with any of several turbulence models. The gradient 
version of this code, which was used for aerodynamic 
sensitivity analysis, was generated by an unconven- 
tional application 31 of the automatic differentiation 
code ADIFOR 32 - 33 to produce a relatively efficient, 
direct mode, gradient analysis code. OFL3D.ADII. 34 
It should be pointed out that the ADIFOR process 
produces a discretized derivative code consistent with 
the discretized function analysis code. The addition of 
a stopping criterion based on the norm of the residual 
of the field equations was the only modification of the 
CFL3D.ADII code made to accommodate the SASDO 
procedure. 
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CFD mesh 
C-O topology 
73x25x25 volume 
49x25 on the wing 



FEM mesh 

3251 elements: 
2240 CST 
1011 truss 

583 points 




Fig. 3 Computational meshes for wing analysis. 

The surface geometry was generated based on the 
parameters described in a previous section by a code 
utilizing the Rapid Aircraft Parameterization Input 
Design (RAPID) technique developed by Smith, et 
al. 35 This code was preprocessed with ADIFOR to gen- 
erate a code capable of producing sensitivity deriva- 
tives as well. 

The CFD volume mesh needed by the flow analysis 
code was generated using a version of the CSCMDO 30 
grid generation code. The associated grid sensitivity 
derivatives needed by the flow sensitivity analysis were 
generated with an automatically differentiated ver- 
sion of CSOMDO. 3 ' In addition to the parameterized 
surface mesh and accompanying gradients. CSOMDO 
requires a baseline volume mesh of similar shape and 
identical topology. The 45,000 grid point baseline vol- 
ume mesh of C-0 topology used in the present flexible 
wing examples was obtained with the Gridgen™ code. 
The wing surface portion of the mesh is shown in 
Fig. 3. This mesh is admittedly particularly coarse 
by current CFD analysis standards. 

The structural analysis code 38 used to compute 
the deflection of the elastic wing was a generic fi- 
nite element code. The flexible structure for the wing 
shown in Fig. 3 was discretized by 583 nodes; there 
were 2,141 constant-strain triangle (CST) elements 
and 1.110 truss elements. Zone boundaries for the de- 
sign variables controlling element size are also shown in 
Fig. 3. Because the elastic deformation was assumed 
to be small, linear elasticity was deemed to be ap- 
propriate. The structural sensitivity equations were 
derived based upon the direct differentiation method. 
We note that the sensitivity of the aerodynamic forces 
appears as a term on the right-hand side (RHS) of the 
deflection sensitivity equations. The derivative of the 
stiffness matrix in these sensitivity equations was also 
generated 31 ' by using the ADIFOR 32,33 technique. We 
note that the coefficient matrix of the structural sensi- 
tivity equations was identical to that of the structural 
equations. Consequently, these structural sensitivity 
equations were solved efficiently by backward substi- 
tution with different RHSs for each sensitivity. 

At the wing surface, i.e., the interface where aero- 
dynamic load and structural deflection information is 
transferred, it was assumed that surface nodes of the 



FEM structural model were a subset of the CFD aero- 
dynamic surface mesh points (see Fig. 3) for the 
present SASDO application. This lack of generality- 
allowed for simplifications in the data transfers and, al- 
though an important issue, it was not deemed crucial 
for these initial SASDO demonstrations. Future ap- 
plications to more complex configurations should allow 
for transfer of conserved information between arbitrary 
meshes as required by the individual disciplines. A 
recent review of such data transfer techniques and a 
specific proposed one are given in Ref. 40. 

Conventional (NAND) and SASDO (SAND) proce- 
dures were implemented using the SQP method of the 
DOT' 11 optimization software. All computations were 
executed on an SCI Origin 2000™ workstation with 
250MHz R10000™ processors. The CFD sensitiv- 
ity calculations were partitioned and run on several 
processors to reduce required memory and elapsed op- 
timization time. This partitioning, however, results in 
additional accumulated computation time due to the 
nature of ADIFOR-generated sensitivity analysis code. 

Results 

The optimization results shown in this work are 
for design problems involving only four or eight de- 
sign variables out of the 21 available wing parameters. 
The results shown by these authors in an earlier work 1 ' 
used design variables that directly affected either the 
aerodynamic analysis alone or both the aerodynamic 
analysis and the structural analysis. In this work, ad- 
ditional design variables are chosen that directly affect 
only the structural analysis. The flow conditions for 
the wing optimizations were M^ = 0.8 and o = 1°. 




Fig. 4 Comparison of planform shapes and sur- 
face pressure contours for 4-design-variable cases, 

M x =0.8. a = 1°. 

Four-Design- Variable Problems 

Table 1 and Fig. 4 show results of several optimiza- 
tion problems involving four design variables: the tip 
chord c t . the tip setback x t , and the structural element 
size factor for the two most inboard zones, F] and 
IV Two of the cases, designated Convc and SASDOr 
in Table 1, represent direct comparisons of SASDO 
and the conventional method for consistent accuracy 
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Table 1 Summary of Four-Design 


-Variable Cases. 




baseline 




Optimizations 








SASDOa 


Convj) 


SASDOc 


Convc 


ft 


1. 


1.07778 


1.07281 


1.07827 


1 .07697 


X( 


1. 


1.93259 


1.91880 


1 ,93052 


1.92614 


n 


1. 


1.06631 


1.08008 


1.09637 


1.09982 


r 2 


1. 


0.69510 


0.69187 


0.68664 


0.68691 


F 


-7.149 


-10.187 


-10.181 


-10.186 


-10.183 


gi 


-.0302 


-.000760 


-.00177 


-.0000629 


-.000915 


g:- 


-.8882 


-0.880 


-0.876 


-0.880 


-0.879 


g3 


-.2647 


-.000639 


-.000159 


-.000149 


+.000108 


C 


9.6e-7 


9.9e-5 


8.9e-7 


6.6e-(5 


4.0e-6 


^ state 


1 


1.1 


4.8 


2.0 


3.9 


Es rad 


— 


27.1 


63.1 


26.9 


40.4 



of function and gradient analyses. The other two cases 
show effects of changes in the accuracy of the function 
and gradient analyses. The resulting designs are es- 
sentially identical for all four cases. Figure 4 shows 
a comparison of the wing planform and the surface 
pressure coefficient results for the baseline design and 
the design designated SASDOc in Table 1 . The shock 
wave has been weakened somewhat in the optimized 
cases from that on the original wing, as one would 
expect. As one can see qualitatively in Fig. 4 and 
numerically from the values of the objective function 
/*", the constraints </,-, and the final design variables 
in fable 1, the final designs are very similar for the 
four problems. The relative computational costs of 



50 r 



m 

u 

o 
!> 

■H 
4J 

id 

H 







conven- 
tional 



SASDO 



Fig. 5 Comparison of computation cost of four- 
design-variable optimization problem using the 
conventional and SASDO methods. 

the optimizations are shown in Table 1 and Fig. 5. 
The accumulated function and gradient analysis times 
(denoted as S.srarf and Ylgrad) are shown separately. 
They have been normalized by the cost of the base- 
line coupled function analysis. The components of the 
analyses were lumped together in Table 1 but they 



have been separated in Fig. 5. The total time for 
performing this optimization problem was reduced by 
36 percent using the SASDO method. The analysis 
alone was reduced by 55 percent, but the gradient eval- 
uation was the dominant cost . 




Fig. 6 Comparison of planform shapes and sur- 
face pressure contours for 8-design-variable cases, 

A/, x = 0.8, a ■ = 1°. 



Eight-Design- Variable Problems 

Table 2 and Fig. 6 show results of three optimization 
problems involving eight design variables: the same set 
used in the four-design-variable eases with the inclu- 
sion of the span 6, the root section max camber z r , 
and the structural element size factor for two more 
zones, F-s and T^. Two of the cases, designated Convjj 
and SASDOb in Table 2, represent direct comparisons 
of SASDO and the conventional method for consis- 
tent accuracy of the function and gradient analyses. 
Figure 6 shows a comparison of wing planform and 
surface pressure coefficient results for the baseline de- 
sign and the design designated SASDOb in Table 2. 
The relative computational costs of the optimizations 
are shown in Table 2 and Fig. 7. The total time for 
performing this optimization problem was reduced by 
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Table 2 Summary of Eight-Design- Variable Cases. 





baseline 


Optimizations 








Conv'A 


SASDOb 


Convy 


Ct 


1. 


0.76574 


0.75929 


0.75794 


Xt 


1. 


2.60000 


2.60000 


2.60000 


b 


1. 


1.05992 


1.07543 


1.07269 


z r 


1. 


1.11511 


1.11173 


1.11175 


Ti 


1. 


2.42444 


2.55263 


2.44498 


r\ 


1. 


2.09537 


2.09924 


2.09394 


r 3 


1. 


1.03635 


1.01883 


1.05778 


r 4 


1. 


0.65137 


0.81046 


0.76612 


F 


-7.149 


-21.255 


-21.235 


-21.244 


gi 


-.0302 


-.7536 


-.7255 


-.7316 


g'-> 


-.8882 


-.00336 


-.00671 


-.000241 


g3 


-.2647 


-.000420 


-.000962 


-.000134 


c- 


9.6e-7 


6.8e-7 


1.8e-5 


l.6e-5 


£ state 


1 


13.5 


4.4 


11.3 


£ s rad 


- 


146.7 


1 15.9 


148.6 



150- 



V 
u 


u 

9 
> 
■H 
V 

id 

H 


m 



100 




n 



d(mesh) 

d(FEM) 

d(CFD) 

mesh 

FEM 

CFD 



conven- SASDO 
tional 

Fig. 7 Comparison of computation cost of eight- 
design-variable optimization problem using the 
conventional and SASDO methods. 

26 percent using the SASDO method. The analysis 
alone was reduced by 60 percent, but, as with the four- 
design- variable problem, the gradient evaluation was 
the dominant cost. 

Further Discussion 

The relative costs, based on CPC timing ratios, 
for the SASDO (SAND) procedures applied to these 
present small 3-D aerodynamic/structural design opti- 



mization problems are about seven-tenths of the costs 
of the corresponding conventional (NAND)procedures. 
This range is very similar to that reported for 2-D 
nonlinear aerodynamic shape design optimization in 
Refs. 1 and 4, even though many of the computational 
details differ. The results given in Ref. 1 were for a 
turbulent transonic flow with shock waves computed 
using a Navier-Stokes code; a direct differentiation ap- 
proach (using ADIFOR) was used for the sensitivity 
analysis. The results reported in Ref. 4 were for a com- 
pressible flow without shock waves computed using a 
nonlinear potential flow code; an adjoint approach was 
used for the sensitivity analysis. Since these two opti- 
mization problems were also not the same, no timing 
comparison between these adjoint and direct differ- 
entiation solution approaches would be meaningful. 
As indicated earlier, an expected speed-up for using 
an adjoint approach instead of the direct differentia- 
tion approach was estimated in Ref. 1. Chattas and 
Bark 1 ' recently reported 2-D and 3-D results for op- 
timal control of steady incompressible Navier-Stokes 
flow that demonstrate an order-of-magnitude reduc- 
tion of CPU time for a SAND approach versus a 
NAN I) approach. These results were obtained using 
reduced Hessian SQP methods that avoid converging 
the flow equations at each optimization iteration. The 
relationship of these methods with respect to other o|>- 
timization techniques is also discussed in Ref. 26. The 
"Control Theory" approach of Jameson" 1 ' and several 
other SAND-like methods for simultaneous analysis 
and design, which were summarized and discussed by 
Ta'asan. 4 - have been applied to aerodynamic shape 
design problems at several fidelities of CFD approx- 
imation. These techniques have obtained an aerody- 
namic design in the equivalent of several analysis CPl 
times for some sample problems. 

Concluding Remarks 

This study has introduced an implementation of the 
SASDO technique for a simple, isolated wing. Initial 
results indicate that SASDO 

1. is feasible under dual simultaneity (i.e. simultane- 
ity not only with respect to analysis and design 
optimization, but also simultaneity with respect 
to flexible wing aero-structural interaction) 

2. finds the same local minimum as a conventional 
technique 

3. is computationally more efficient than a conven- 
tional gradient-based optimization technique 

4. requires few modifications to the analysis and sen- 
sitivity analysis codes involved 

5. is effective at reducing the function analysis cost, 
but the gradient analysis time is the dominant 
cost 



9 OF 11 



American Institute of Aeronautics and Astronautics Paper 2001-2527 



Acknowledgment 

The second author, G. J.-VV. IIou, was supported in 
this work by NASA through several Tasks under con- 
tract NASI- 19858 and NASA P.O. No. L-9291 with 
the Old Dominion University Research Foundation. 

References 

'Newman. P. A.. Hou. G. J.-W.. and Taylor, III, A. C, 
"Observations Regarding Use of Advanced CFD Analysis, Sen- 
sitivity Analysis, and Design Codes in MIX)," Multidisciplinary 
Design Optimization: State of tht Art. edited by N. M. Alexan- 
drov and M. Y. Hussaini. SIAM Proceedings Series, Philadel- 
phia, 1997. pp. 263 279. also IOASE Report 96-16, NASA CR 
198293, (available electronically at www.icase.edu). 

2 Newman, III, J. C, Taylor. Ill, A. C, Barnwell, R. W., 
Newman, P. A., and Hou, G. J.-W.. "Overview of Sensitivity 
Analysis and Shape Optimization for Complex Aerodynamic 
Configurations." Journal of Aircraft. Vol. 36, No. 1, 1999. 
pp. 87-96. 

Alexandrov, N. M. and Hussaini, M. Y.. editors. Mul- 
tidisciplinary Design Optimization: State of the Art, SIAM 
Proceedings Series, Philadelphia. 1997. 

4 Ghattas, O. and Orozco, C. E., "A Parallel Reduced Hes- 
sian SQP Method for Shape Optimization," Multidisciplinary 
Design Optimization: State of the Art, edited by N. M. Alexan- 
drov and M. Y. Hussaini, SIAM Proceedings Series, Philadel- 
phia. 1997. pp. 133-152. 

Tlou, G. J.-W., Taylor, III, A. C, Mani, S. V., and New- 
man, P. A., "Simultaneous Aerodynamic Analysis and Design 
Optimization," Abstracts from 2nd L'.S. National Congress on 
Computational Mechanics, Aug. 1993, pp. 130. 

'Mani, S. V.. Simultaneous Aerodynamic Analysis and De- 
sign Optimization, Master's thesis. Old Dominion University, 
Norfolk, VA, Dec. 1993. 

7 Hou. G. J.-W., Korivi. V. M.. Taylor, III. A. O. Maroju, 
\'., and Newman. P. A., "Simultaneous Aerodynamic Analysis 
and Design Optimization (SAADO) of a Turbulent Transonic 
Airfoil Using a Navier-Stokes Code With Automatic Differen- 
tiation (ADIFOR)." Proceedings of the Computational Aero- 
scu rices Workshop 95. edited by W. J. Feiereisen and A. K. 
Lacer. NASA CD OP-20010. Jan. 1996. pp. 82-85. 

*Gumbert. C. Ft.. Hou, G. J.-W., and Newman, P. A., 
"Simultaneous Aerodynamic Analysis and Design Optimiza- 
tion (SAADO) of a 3-D Rigid Wing," Proceedings of 14th 
AIAA Computational Fluid Dynamics C'onference, Norfolk, 
June 1999. pp. 402-118, also AIAA Paper 99-3296. 

9 Gumbert. C. R.. Hou. G. J.-W., and Newman, P. A., 
"Simultaneous Aerodynamic Analysis and Design Optimization 
(SAADO) of a.3-D Flexible Wing," AIAA Paper 2001-1 107, Jan. 
2001. 

( Haftka, R. T., Sobieszczanski-Sobieski, J., and Padula, 
S. L.. "On Options for Interdisciplinary Analysis and Design 
Optimization," Structural Optimization. Vol. 4. 1992, pp. 65- 
74. 

"Shubin, G. R., "Application of Alternative Multidis- 
ciplinary Optimization Formulations to a Model Problem 
for Static Aeroelasticity." Journal of Computational Physics, 
Vol. 118, No. 1. 1995. pp. 73 85. 

"Arian, E., "Analysis of the Hessian for Aeroelastic Opti- 
mization," ICASE Report 95-84, Dec. 1995. 

''Walsh, J. L.. Townsend. J. C, Salas, A. O., Samareh. J. A.. 
Mukhopadhyay, V., and Barthelemy. J. F., "Multidisciplinary 
High! idelity Analysis and Optimization of Aerospace Vehicles, 
Part 1: Formulation," AIAA Paper 2000-0418, Jan. 2000. 

14 Walsh, J. L.. Weston, R. P., Samareh, J. A., Mason, B. 
H.. Green. L. L.. and Biedron. R. T.. "Multidisciplinary High- 
Fidelity Analysis and Optimization of Aerospace Vehicles, Part 
2: Preliminary Results." AIAA Paper 2000-0419. Jan. 2000. 



15 Reuther. J. J., Alonso, J. J., Martins, J. R., and Smith, 
S. C, "A Coupled Aero-Structural Optimization Method for 
Complete Aircraft Configurations," AIAA Paper 99-0187. Jan. 
1999. 

16 Kapania, R. K., Eldred, L. B., and Barthelemy. J. -F. M.. 
"Sensitivity Analysis of a Wing Aeroelastic Response," Journal 
of Aircraft, Vol. 30, No. 4, 1993. pp. 496-504, also AIAA Paper 
91-1103. 

''Arslan, A. E. and Carlson, L. A., "Integrated Determi- 
nation of Sensitivity Derivatives for an Aeroelastic Transonic 
Wing," Proceedings for the 5th A1AA/USAF/NASA/ISSM0 
Symposium on Multidisciplinary Analysis and Optimization. 
Panama City Beach, FL, Sept. 1994, pp. 1286-1300, also AIAA 
Paper 94-4400 CP. 

l8 Ciunta, A. A. and Sobieszczanski-Sobieski, J.. "Progress 
Toward Using Sensitivity Derivatives in a High-Fidelity Aero- 
elastic Analysis of a Supersonic Transport," Proceidings of 
the 1th AIAA/VSAF/NASA/1SSMO Symposium on Multidis- 
ciplinary Analysis and Optimization. St. Louis, MO, Sept. 1998, 
pp. 441-453. also AIAA Paper 98-4763 CP. 

19 Giunta, A. A., "Sensitivity Analysis for Coupled Aero- 
Structural Systems." NASA TM- 1999-209367, Aug. 1999. 

20 Newman, III, J. O. Whitfield, D. L„ and Anderson, K. 
W., "A Step-Size Independent Approach for Multidisciplinary 
Sensitivity Analysis and Design Optimization." Procttdings of 
the l?th Applied Aerodynamics C'onference. Norfolk, VA. June 
1999, pp. 12 22. also AIAA Paper 99-3101. 

2l Ghattas, O. and Li. X., "Domain Decomposition Methods 
for Sensitivity Analysis of a Nonlinear Aeroelastic Problem." In- 
ternational Journal of Computational Fluid Dynamics, Vol. 11. 
1998, pp. 113-130. 

2 -M0ller, H. and Lund. E., "Shape Sensitivity Analysis 
of Strongly Coupled Fluid-Structure Interaction Problems." 
Procttdings of 8th AIA A/VSAF/N ASA/ISSMO Symposium 
on Multidisciplinary Analysis and Optimization. Long Btach. 
Sept. 2000, also AIAA Paper 2000-4823. 

2! Maute, K., Nikbay, M.. and Farhat, C, "Analytical Based 
Sensitivity Analysis and Optimization of Nonlinear Aeroelas- 
tic Systems." Proceedings of 8th AIAA/USAF/NASA /ISSMO 
Symposium on Multidisciplinary Analysis and Optimization, 
Long Beach, Sept. 2000. also AIAA Paper 2000-4825. 

24 Hou. G. and Satyanarayana, A., "Analytical Sensitivity 
Analysis of a Static Aeroelastic Wing." Proceedings of 8th 
A1AA/USAF/NASA/ISSMO Symposium on Multidisciplinary 
Analysis and Optimization, Long Beach, Sept. 2000, also AIAA 
Paper 2000-4824. 

25 Biros, G. and Ghattas. O., "Parallel Newton- Krylov Al- 
gorithms For PDE-Constrained Optimization," Proceidings of 
Scientific Computing 99. Portland, OR, 1999. 

26 Ghattas, O. and Bark, J.-H., "Optimal Control of Two- 
and Three-Dimensional Navier-Stokes Flows," Journal of Com- 
putational Physics. Vol. 136. No. 2, 1997, pp. 231 244. 

-''Pierce, N. A. and Giles. M. B., "Adjoint Recovery of Su- 
perconvergent Functionals From PDE Approximations." SIAM 
Review, Vol. 42, No. 2, 2000, pp. 247-264. 

2e Venditti,D. A. and Darmofal.D. L., "Multilevel Error Esti- 
mation and Grid Adaptive Strategy for Improving the Accuracy 
of Integral Outputs," Proceedings of the AIAA 15th Compu- 
tational Fluid Dynamics Conference, Norfolk, June 1999. also 
AIAA Paper 1999-3292. 

29 Jameson, A. and Vassberg. J. O, "Computational Fluid 
Dynamics for Aerodynamic Design: Its Current and Future Im- 
pact." AIAA Paper 2001-0538, Jan. 2001. 

M Rumsey. C Biedron, R., and Thomas. J., "CFL3D: Its 
History and Some Recent Applications," NASA TM-112861, 
May 1997. 

31 Sherman, L., Taylor, III, A., Green, L., Newman, P.. Hou, 
G., and Korivi, M., "First- and Second- Order Aerodynamic 
Sensitivity Derivatives via Automatic Differentiation with lncre- 



1(1 OF 11 



American Institute of Aeronautics and Astronautics Paper 2001-2527 



mental Iterative Methods," Journal of Computational Pliysics. 
Vol. 129, No. 2, 1996, pp. 307 336. 

32 Bischof. C. H., Carle. A., Corliss. G. F.. Griewank, A., 
and Hovland. P., "ADIFOR: Generating Derivative Codes from 
Fortran Programs," Scientific Programming. Vol. 1, No. 1, 1992, 
pp. 1-29. 

"Bischof, C. and Griewank. A., "ADIFOR: A Fortran Sys- 
tem for Portable Automatic Differentiation," Proceedings of tht 
4th AIAA/USAF/NASA/OAI Symposium on Multidisciplinary 
Analysis and Optimization. Cleveland, OH, Sept. 1992. pp. 433- 
441. also AIAA Paper 92-4744 CP. 

34 Taylor, III. A. C, Oloso, A., and Newman, III, J.C., 
"CFL3D.ADII (Version 2.0): An Efficient. Accurate. General- 
Purpose Code for Plow Shape-Sensitivity Analysis." AIAA Pa- 
per 97-2204, June 1997. 

3 ' Smith, R. E., Bloor. M. I. G., Wilson, M. J., and Thomas, 
A. T.. "Rapid Airplane Parametric Input Design (RAPID)," 
Proceedings of the 12th AIAA (Computational Fluid Dynamics 
Conference . San Diego, June 1995. pp. 452-462. also AIAA Pa- 
per 95-1687. 

36 Jones. VV. T. and Samareh-Abolhassani, J., "A Grid Gener- 
ation System for Multidisciplinary Design Optimization." Pro- 
ceedings of tht 12th AIAA Computational Fluid Dynamics Con- 
ferenct. San Diego. June 1995, pp. 474-482. also AIAA Paper 
95-1689. 

37 Bischof, C Jones. \V. T., Samareh-Abolhassani. J., and 
Mauer, A., "Experiences with the Application of the AD1C Au- 
tomatic Differentiation Tool to the CSCMDO 3-D Volume Grid 
Generation Code," AIAA Paper 96-0716. Jan. 1996. 

38 Nguyen. D. T., "Finite Element Software for Multidisci- 
plinary Design Optimization. Final Report," Tech. rep.. Old Do- 
minion University Research Foundation. Contract NASI- 19858, 
Task 69. NASA Langley Research Center. Hampton. VA. Nov. 
1995. 

39 Hou. G., Arunkumar, S.. and Tiwari. N. S., "First- and 
Second-Order Sensitivity Analysis of Finite Element Equa- 
tions via Automatic Differentiation," Proceedings of tht 7th 
AIAA/USAF/NASA/1SSMO Symposium on Multidisciplinary 
Analysis and Optimization, St. Louis, MO. Sept. 1998, pp. 454- 
464. also AIAA Paper 98-4764. 

40 Samareh, J. A. and Bhatia, K. G.. "A Unified Approach to 
Modeling Multidisciplinary Interactions." Proceedings of the 8th 
AI A A/NASA/US A F/ISSMO Symposium on Multidisciplinary 
Analysis and Optimization. Long Beach, Sept. 2000, also AIAA 
Paper 2000-4704. 

'"Anon.. DOT ('sers Manual: Version 5.0. Vanderplaats 
Research & Development. Inc.. Colorado Springs, .May 1999. 

42 Ta'asan, S., "Trends in Aerodynamic Design and Optimiza- 
tion: A Mathematical Viewpoint." Proceedings." Proceedings of 
tht 12th .47.4.4 Computational Fluid Dynamics Conference . San 
Diego. June 1995. pp. 961-970. also AIAA Paper 95-0001. 



11 OF 1 1 

American Institute ok Aeronautics and Astronautics Paper 2001 2527 



